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PACS 73 . 50 . Lw - Thermoelectric effects 
PACS 05 . 60 . Gg - Quantum transport 

PACS 73.63.Kv - Electronic transport in nanoscale materials and structures: QD 

Abstract -We study linear response and nonequilibrium steady-state thermoelectric transport 
through a single-level quantum dot tunnel coupled to two reservoirs held at different temperatures 
as well as chemical potentials. A fermion occupying the dot interacts with those in the reservoirs 
by a short-ranged two-particle interaction. For parameters for which particles flow against a 
bias voltage from the hot to the cold reservoir this setup acts as an energy-conversion device 
with which electrical energy is gained out of waste heat. We investigate how correlations affect its 
efficiency and output power. In linear response the changes in the thermoelectric properties can be 
traced back to the interaction induced renormalization of the resonance line shape. In particular, 
small to intermediate repulsive interactions reduce the maximum efficiency. In nonequilibrium the 
situation is more complex and we identify a parameter regime in which for a fixed lower bound of 
the output power the efficiency increases. 



Introduction. — Quantum dots (QDs) arc considered 
as high-potential solid-state energy conversion devices to 
gain electrical energy out of waste heat |T|. In a minimal 
setup electrons from a hot right reservoir (at temperature 
' Tr) are transported via a weakly coupled QD to a cold left 
[ one (at T^) against a bias voltage V ^ ji^ — fijf > 0. Here 
I fia denotes the chemical potential of reservoir a = L, R 
] and we use units in which the elementary charge e = 1. 
• Within linear response, that is for small temperature dif- 
I ferenccs and small bias voltages, one aims at a high effi- 
' ciency rj of energy conversion. It is defined as 77 = P/|/h|, 
'with the output power P = V\Ic\ and the charge/heat 
I currents Ic/h [2] ■ As a measure of rj in QD setups one 
often uses the so-called figure of merit ZT [T]; the larger 
ZT the closer 77 comes to its upper bound, the Carnot effi- 
ciency 77c . We here neglect the contribution of the phonon 
thermal conductance to ZT and focus on the electronic 
degrees of freedom. The figure of merit is inversely pro- 
portional to the Lorenz number. Under the assumption of 
scattering processes which on the scale Tr sa w T only 
weakly depend on energy the Wiedemann- Franz law states 
that the latter is a material independent number. This in 
turn limits the electronic ZT in bulk materials which obey 
this law. In QDs due to the strong energy dependence of 
the transmission close to transport resonances [T1[31I3] the 
Wiedemann-Franz law does not apply as long as the tem- 



perature is larger than the resonance width F. Therefore, 
in transport through dots with sharp resonances a large 
figure of merit and thus a high efficiency can be realized 
[UlSHSI. Equivalently, QDs can be used as efficient solid- 
state systems for cooling of nanoelectronic devices [6]. 

Due to the confinement of electrons in QDs to meso- 
or nanoscopic regions the local Coulomb interaction be- 
comes a relevant energy scale. This two-particle interac- 
tion strongly alters the line shape of resonances and thus 
the thermoelectric properties [4l[5l[7]. For example, the 
Kondo effect, resulting from correlated spin fluctuations, 
leads to a very sharp many-body resonance of width much 
smaller than the noninteracting F. This has crucial con- 
sequences for thermoelectric transport [8Hl4] . Even in the 
absence of the Kondo effect correlated charge fluctuations 
affect the line width as was worked out in detail for the in- 
teracting resonant level model (IRLM) [15H20] . For weak 
to intermediate local Coulomb repulsions the interaction 
renormalized width F''*'" is increased compared to the non- 
interacting one F, while for strong repulsive interactions as 
well as attractive ones F''°" < F. Effectively attractive in- 
teractions might be realized in molecular QDs with strong 
local electron-phonon coupling The consequences of 
the interaction renormalized resonance line shape for the 
thermoelectric transport properties were not discussed so 
far. We fill this gap and show that the parameter de- 
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pendence of 77 and P can be understood in terms of the 
line width renormalization. We provide approximate an- 
alytical results for the efficiency and the output power 
which for small two-particle interactions agree very well 
with numerical data. For weak repulsive interactions the 
maximum efficiency decreases, it increases for attractive 
ones. Depending on the temperature regime considered 
the two-particle interaction can lead to an increase or to 
a decrease of the power at maximum efficiency. 

In linear response the output power is of order AT^, 
with AT = Tl - Tr < 0, and thus small (see below). 
Depending on the precise technological conditions under 
which the QD 'heat engine' is supposed to perform (fixed 
or tunable AT, limited or unlimited supply of heat, scal- 
ability in parallel or series, etc.) it might be meaningful 
to maximize P instead of r] by varying the model param- 
eters and the voltage at fixed temperatures, with |AT| 
not necessarily being small. In fact, within the frame- 
work of nonequilibrium thermodynamics studying the ef- 
ficiency at maximum power of general heat engines is an 
active field of current research [2TH27] . In this context the 
noninteracting resonant level model in the limit of van- 
ishing level-reservoir coupling was studied as a toy model 
|28| . Even questions such as 'what is the optimal efficiency 
reachable for the power being larger than a given lower 
bound' might be of interest. Computing the efficiency 
at maximum power and answering questions of the above 
type requires access to the full nonequilibrium (in V and 
AT) steady-state properties |29| . For microscopic models 
of QDs with local two-particle correlations achieving this 
presents a formidable challenge; for a recent review, see 
e.g. Ref. |30| . In particular, this holds if nonperturbative 
effects in both the level- reservoir coupling as well as the 
two-particle interaction become crucial as it is the case in 
the IRLM [T5H20]. 

Only recently a very flexible nonperturbative tool was 
developed which allows to treat the full-fledged nonequi- 
librium steady state of the IRLM [T9ll3ip32j at weak inter- 
actions. It is based on Keldysh Green functions and 
the functional renormalization group approach to quan- 
tum many-body physics |34| . Here we apply this method 
to study the efficiency and power output of the IRLM 
'heat engine' beyond linear response. Broadly speaking 
the effect of the interaction on the efficiency (at maxi- 
mum power) is similar to the one found in linear response 
(summarized in the next to last paragraph) . By closer in- 
spection we identify a situation in which a weak repulsive 
interaction can lead to an increase of the efficiency by a 
few percent at fairly large power output. 

We here focus on correlation effects in the elementary 
two-reservoir IRLM but expect our results to be of rele- 
vance for other models showing correlated charge fiuctu- 
ations as well. Very recently three-terminal systems with 
the reservoirs having more complex degrees of freedom 
(e.g. magnetic ones) were identified as promising energy 
converters [35H37] . In those systems electrons are trans- 
ported against a bias between two of the reservoirs by 
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Fig. 1: Sketch of the investigated model. 
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extracting heat from the third one. 

Quantum dot modeL — Our model Hamiltonian 
sketched in Fig.[T]consists of three parts H = H^+Hd+Hc- 
The two spinless reservoirs are assumed to be identical and 
are given by 

= Xl^'=4,aCfe,a> (1) 

with fermionic ladder operators c^^^ , where k denotes a set 
of quantum numbers characterizing the reservoir states. 
To disentangle effects of the reservoir band structure en- 
coded in the dispersion tk and the correlations on the 
thermoelectric properties we assume structureless reser- 
voirs with an energy independent density of states p be- 
tween ~D and D and the bandwidth 2D being much larger 
than any other energy scale. The reservoirs arc in grand 
canonical equilibrium with T^ = l//3a and chemical po- 
tentials centered aroimd zero: /.j^ = — Mfl — V/2 (with 
h = = e = I). The dot part of H is given by ~ en, 
with the dot occupation number operator n = Sd. The 
energy e of the level is assumed to be tunable; in experi- 
ments this is achieved by applying a voltage to a properly 
fabricated gate. Finally, the coupling reads 
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where : . . . : denotes normal ordering. This and the shift 
of the dot occupancy by —1/2 ensures that e = corre- 
sponds to half dot filling. For simplicity we assume the 
tunnel couplings to the left and right reservoir t to be 
equal; similarly for the two-particle interactions u to the 
left and right. We note that within our approach to the 
many-body problem these restrictions can easily be re- 
laxed. We use the dimensionless interaction U = pu. 

Thermoelectric transport at [/ = 0. For U = 

transport through the dot is characterized by a Breit- 
Wigner transmission resonance 



{lo - e)2 + r2 



(3) 



of width r = 2-11 pt^. The steady-state charge current Ic 
(leaving the left reservoir) and heat current /h (entering 
the right reservoir) can be computed using the Landauer- 
Biittiker formalism 



(4) 



I duj {uj- ij,r)t{uj) [/l(w) - /i?(w)] , (5) 
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Fig. 2: (Color online) Output power P of our quantum dot 
device for vanishing two-particle interaction as a function of 
e/r and V/T for Tl/T = 1 and Tr/T = 20. Data points are 
only shown in the regime in which the device acts as a 'heat 
engine' (curved surface with lines). 



with the Fermi functions /^(a;) = [e('^-'^°)/^<» l] .We 
express the energy integral in /c in terms of the digamma 
function ^'(0, z) as 



/e(e) = ;^Im 

ZTT 



(6) 



and the heat current as a sum of the charge current and 
the Hilbert transform of the latter 



^h(e) = (e - fJ-R)Ic{e) - 



de 



(7) 



In Eqs. ^ and ([7|) F can be scaled out (taken as the unit 
of energy) and rj as well as P can be investigated as func- 
tions of the model parameter e/F, the external bias voltage 
V/T, and T^i^lV. In analogy to studies on periodic heat 
engines in thermodynamics we assume that T^/^j are fixed 
by the environment. With |Ar|/F not necessarily being 
small we cannot resile to the linear response regime. Fig- 
ured] shows P = T^|/c| in the e-T^-plane for Tl/F = 1 and 
Tr/F = 20. It has a unique maximum. 

To set the stage for more general considerations we first 
consider the limit F — )- The second term in Eq. ([7]) 
is of higher order in F as compared to the first one and 
can be neglected. As a consequence charge and heat flow 
are 'perfectly coupled' (proportional to each other); for 
F — > the energy of every particle coming from the dot 
level becomes sharp and is fixed at e [55] . For applications 
it is meaningful to maximize the output power — the max- 
imum is denoted by Pm in the following — and study the 
efficiency at maximum power t]^^ as a function of Tj^/n. 
The optimization is performed with respect to the level 
energy e and the externally applied voltage V both being 
parameters which in experiments on QDs can routinely 
be varied with high precision. As shown analytically in 
Ref. [IH] foi' F 0, 77nip and Pm are functions of Tl/Tr 
only and can thus be written as functions of the Carnot 
efficiency r?c = 1 - Tl/Tr = \AT\/Tr. 




'/mp 




0.7 






0.6 




0.8 


0.5 






0.4 




0.6 


0.3 




O 

55-0.4 


0.2 
0.1 


■ 

■ 


0.2 
0.0 




Fig. 3: (Color online) The efficiency at maximum power (left) 
and the maximum power (right) for U = Q as. functions of Tl /F 
and Tjc = 1 ~ Tl/Tr. Contour lines are shown. 



We now return to F > and numerically maximize the 
output power with respect to e and V using Eqs. ^ and 
([7]). Both, 77mp and P,„ are no longer functions of rjc only. 
This becomes apparent from Fig. [3] which shows r/,„p and 
Pni as functions of Tl /F and 77c ■ Only in the limit of large 
Tl/T, that is small F, the dependence on Tl/T drops out. 
For small Tl/T rather large rye = |AT|/Tr (close to 1) are 
required to obtain a sizable efficiency at maximum power. 

Expanding Eqs. ([6]) and ([7]) to lowest order in AT/T and 
V/T we obtain linear response results [5] . The currents fol- 
low from the Onsager matrix containing the conductances 



Ic 

h 



Gh 



Gl 



V 

AT 



(8) 



The Onsager (time-reversal) symmetry gives Gj^ — 
TG^"^ , with T = Tl = Tr, and the independent matrix 
elements can be written as 



Gc 



* 1 



TT + W 



G: 



AT 



Gh 



iT 
^Stt^T 
T /T 



1, 



2tt 

t: + w 



27r 



2tt \T 



1 

An 



w^-^ 1, 



IT + W 

2tt 



(9) 

(10) 

(11) 



with the trigamma function ^'(1, z) and w = (F + ie)/T . 

In linear response one aims at a high efficiency. In- 
dependent of the model considered on first maximizes 
•q = V\Ic\/\h\ with respect to F [22]25] using Eq. dU 
and leading to 



^ZT +1-1 lATl 



VZT +1+1 T 



ZT 



GcGh-(G^^)2T- 



(12) 

The maximum is reached for a voltage of order AT and 
the linear response regime is not left. Equation (|12p con- 
stitutes the relation between the efficiency and the figure 
of merit ZT announced in the introduction; the latter be- 
ing expressible in terms of the conductances. We can now 
be more precise: ZT is a measure for the linear response 
efficiency maximized with respect to V . For ZT — > 00, 770 
approaches the linear response Carnot efficiency \AT\/T 
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Fig. 4: (Color online) Linear response maximum efficiency (up- 
per panel) and the corresponding output power (lower panel) 
for different interactions U as functions of T/F. The lines fol- 
low from Eqs. (|8}-(IlT} and (for C/ / 0) while the symbols 
are numerical data obtained from our many-body approach. 
Inset: zoom-in of the regime T/F ^ 1. 



which constitutes an upper bound for all thermodynamic 
efficiencies. The power corresponding to r/o is given by 



Pn = 



TGh vzr+T- 1 /AT 



^/ZT+l ^JZT +1 + 1 V 



(13) 



and thus of order (AT)^; one factor AT comes from the 
voltage at which the maximum efficiency is reached and 
the other from the charge current. In the limit ZT — > oo, 
with rjo ^ r]c, the prefactor of Pq vanishes. 

We now return to our specific model and similar to the 
nonequilibrium case above first consider the limit T ^ 0. 
Independent of e ^ one finds ZT ^ l/F [5] and thus 
Vo Vc U&i- Furthermore, Gh ^ F and Pq 0. 

For F > we numerically maximize tjq (in addition) 
with respect to e using Eqs. ([51)- PT|) [5]. The resulting 
efficiency rj^ and power at maximum efficiency Pmc as 
functions of T/F are shown in Fig. As discussed above 
for F — )■ 0, that is for large T/F, r^m approaches Carnot 
efficiency and Pmc/AT^ vanishes. The latter also vanishes 
for T/F — >■ and is thus a nonmonotonic function with a 
maximum at T sa F. Our calculations show that for small 
T the level energy at which optimal efficiency is reached 
approaches a constant. Using this and the expansions of 
Eqs. (©-ini) wc obtain TGh and ZT - T^. Thus 



P„ic/AT2 - T2 (see the inset of Fig. 13). 

Thermoelectric transport at C/ 7^ 0. To investi- 
gate the linear response and nonequilibrium steady-state 
transport properties in the presence of two-particle cor- 
relations we use a technique which is based on nonequi- 
librium Keldysh Green functions |33| and the functional 
renormalization group approach to quantum many-body 
physics j34j . In this one derives an exact infinite hierar- 
chy of coupled differential equations for the many-body 
self-energy and the effective n-particle interactions (one- 
particle irreducible vertex functions). The derivative is 
taken with respect to a low-energy cutoff and the set 
of equations have to be integrated from infinity to zero. 
Concrete calculations require truncations. We apply the 
lowest-order approximation in which only the self-energy 
flow is kept. It was earlier used to study the IRLM for 
\U\ <C 1 in equilibrium as well as the charge current in 
steady-state nonequilibrium [19l[3Tl[32] . We refrain from 
giving any technical details and instead refer the interested 
reader to Refs. [TH] and |32j . 

As its main effect the two-particle interaction gener- 
ates a renormalization group flow of the hybridization 
F/2 to the individual reservoirs. For the left one it is 
cut off by the largest of the energy scales |e — Tl, 
and the total hybridization; for the right one one has to 
replace the index T by P [HI [311132]. This leads to a 
power-law scaling F™"/F ^ {s/D)~'^^'^\ with the corre- 
sponding largest scale s (still much smaller than D) and 
the f7-dependent exponent v. To leading order one finds 
iy(U) = 2U ITil[IHHkl[3Ill31[3i- This correlation effect 
cannot be captured by perturbative (in either U or F) ap- 
proaches. In a general nonequilibrium setup F^f" 7^ FJf" 
even in the case of equal bare hybridizations. To exemplify 
how the renormalization manifests in physical observables 
we consider the linear response (F, AT — )• 0) charge con- 
ductance Gc as a function of e close to the resonance (at 
e = 0) and T = [T71[T1[1D]. The largest scale cutting 
off the renormalization group flow of the hybridizations is 
given by F . Using the above scaling law one concludes 
that the resonance in Gc has the width 



2F 



-2U 



(14) 



which deflnes the emergent energy scale W. The reso- 
nance becomes wider for weak (to intermediate [17]) re- 
pulsive interactions but more narrow for attractive ones. 
To illustrate how F''°" changes when two or more of the 
energy scales are of the same order [TOIIHTIIS^ we consider 
T >W. In this case F"-™ « W(T/W)-^^; see Fig. 1 of 
Ref. [32]. Since |J7| ^ 1 the correction to W is of order 
one as long as T does not exceed W by several orders of 
magnitude. Further down this will become crucial. 

We start the discussion of our results of correlation ef- 
fects on the efficiency and power considering linear re- 
sponse. The symbols in Fig. |4] show the maximum ef- 
ficiency T]rn (maximized with respect to V and e) and 



p-4 



EfRciency and power of a thermoelectric quantum dot device 



the power at maximum efficiency Pmc for U = ±0.1/7r 
obtained by numerically integrating the renormalization 
group equations and computing the charge and thermal 
currents |39j . For weak repulsive interactions -q^^ is re- 
duced compared to the noninteracting result, while it is 
increased for weak attractive ones. Depending on the tem- 
perature regime considered the interaction enhances or re- 
duces Pme- These results can be understood analytically 
by employing the above discussed energy scale dependence 
renormalization of the resonance width. For [7 = the 
level energy e at which the maximum efficiency is reached 
is of the order of a few T at large T [S] while it saturates 
at a constant of order F for small T . This also holds for 
small \ U\. Thus Eq. p4)) gives a very good estimate of the 
renormalized hybridization for all temperatures shown in 
Fig. 21 This motivates us to replace F in the noninter- 
acting expressions Eqs. ([5])- PT|) for the charge and heat 
currents by W Eq. ([Ti)) . The corresponding results for 
and Pme are shown as dashed lines in Fig. |4l As expected 
they show excellent agreement with the numerical data. 
The deviations at T/F ^ 1 are consistent with the finite 
T correction to W discussed above. 

As our approximate approach is restricted to \U\ <C 1 
we cannot access the regime of large interactions. It is 
known that F''°" is a nonmonotonic function of ?7 > and 
for large repulsive interactions one finds F''°" < F [TT l lTSj. 
This can be explained by the appearance of a second order 
term ^ —U^ (related to the so-called orthogonality catas- 
trophe) which dominates for [/ > 1 and changes the sign 
of the exponent v [18]. Based on our above results it is 
thus reasonable to assume that for large repulsive interac- 
tions 77 will be increased compared to the noninteracting 
case. It would be very interesting to explicitely verify this 
using a method complementary to ours. 

We next leave the linear response regime. Figure [5] 
shows the U dependence of the maximum power Pm and 
the efficiency at maximum power Tynip for different sets 
of Tl/_r obtained numerically using our nonequilibrium 
many-body method. In analogy to the maximum effi- 
ciency in linear response Tymp is a decreasing function of 
U . A similar analogy to the linear response results for Pme 
Fig. |4] holds for the maximum power Pm- It is a decreasing 
function of U for temperatures Tl/_r smaller than F and 
an increasing one for Ti^j^jV ^ 1. For intermediate tem- 
peratures we find nonmonotonic behavior. We have also 
studied the JJ dependences of the maximum (with respect 
to V and e) efficiency and the corresponding output power 
beyond the linear response regime. They are qualitatively 
similar to the results for T^mp and Pm shown in Fig. [5] 

We finally investigate a situation in which even a weak 
repulsive interaction can lead to an increase of the effi- 
ciency. Let us assume that our 'heat engine' has to pro- 
duce at least a certain minimal output power. This might 
e.g. be a sensible requirement if we are interested in charg- 
ing a battery. For appropriate temperatures Pi/fl the 
maximum power increases with V (see Fig. [5]) and one 
can hope to overcome the decrease of the efficiency. For 
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Fig. 5: (Color online) Efficiency at maximum power (upper 
panel) and the maximum power (lower panel) for different sets 
of r^/fl as functions of U . The lines are guide to the eyes. The 
data are scaled by the given factors. 



fixed Tj, jji we proceed as follows: as our lower bound of the 
power Pbound we take the maximum power at J7 = 0. This 
is a convenient choice to reduce the number of parameters. 
Varying e and V we then search for the largest efficiency r\y 
at fixed [/ > with P > Pbound- As this is a numerically 
demanding procedure — remind that for every parameter 
set the currents /c/h have to be computed numerically 
solving coupled differential flow equations — we restrict 
ourselves to one of the temperature sets (with increasing 
power) of Fig. El namely Tl/F = 1 and Tr/F = 20. The U 
dependence of ry> is shown in Fig. |6l Under the above re- 
quirement one can increase the efficiency by a few percent 
by turning on a small repulsive interaction. 

Summary. — Wc have discussed the charge and heat 
transport properties of an elementary single-level quan- 
tum dot 'heat engine' which can be used to convert waste 
heat into electrical energy. Firstly, wc have studied the 
linear response and nonequilibrium steady-state thermo- 
electric properties in the absence of local Coulomb correla- 
tions. The efficiency and output power of this device was 
investigated in all details. We next included the local two- 
particle interaction which in meso- or nanoscopic devices 
is an important energy scale. Using a flexible approxi- 
mate many-body method which can be applied in equilib- 
rium (linear response) as well as the in the nonequilibrium 
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Fig. 6: (Color online) Interaction U dependence of the largest 
efficiency obtainable with P > Pbound for Tl/T = 1 and 
Tr/F = 20 (for details see the text). The line is an Akima 
spline interpolating between the computed data points (sym- 
bols). The efficiency can be raised by a few percent by turning 
on a small repulsive interaction. 



steady state we were able to present a comprehensive pic- 
ture of the correlation effects on the efficiency and output 
power. 
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